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Abstract 

Within a transport code simulation for heavy-ion collisions at bombarding 
energies around 1 AGeV, we demonstrate that double-differential di-electron 
spectra with suitable kinematical cuts are useful to isolate (i) the p meson 
peak even in case of strong broadening, and (ii) the in-medium u decay con- 
tribution. The expected in-medium modifications of the vector meson spectral 
densities can thus be probed in this energy range via the di-electron channel. 
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I. INTRODUCTION 

One of the main goals of the experiments with the HADES detector [1] at the heavy- 
ion synchrotron SIS18 at GSI/Darmstadt is the search for a direct evidence of in-medium 
modifications of hadrons via the di-electron (e + e _ ) channel. This goal is supported by 
various theoretical indications concerning an important sensitivity of vector mesons to partial 
restoration of chiral symmetry in a dense and hot nuclear medium [2]. In particular, as shown 
within lattice QCD [3], the chiral quark condensate (qq) as order parameter of the chiral 
symmetry breaking decreases with increasing baryon density and temperature. Since (qq) 
is coupled to the vector meson spectral densities one can expect considerable in-medium 
effects of light vector mesons which manifest themselves, e.g., as resonance mass shifts and 
broadening of widths [4] . 

In spite of the efforts spent to predictions [5-13], there is no commonly accepted and re- 
liably quantitative estimate of the change of the mentioned parameters or spectral functions 
in a nuclear medium. For instance, within the QCD sum rule analysis [5-7] the in-medium 
modification of the vector meson masses in zero-width approximation is dictated by the 
density and temperature dependence of the quark condensates. While the two-quark con- 
densate (qq) can be estimated by a model independent approach with reasonable accuracy, 
the so far poorly known four-quark condensates like ({qq) 2 ) lead to uncertainties for the 
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prediction of the vector meson mass shift [14] . From this point of view, direct measurements 
of the vector meson masses and widths, as being at reach with the high-resolution HADES 
detector system [1], can also shed light on fundamental problems such as the in-medium 
change of the four-quark condensate. 

Since the most pronounced vector meson mass shift is expected for large baryon density 
in the nuclear medium, there is the hope to verify the effect in medium-energy heavy-ion 
collisions at SIS 18 energies, where baryon densities up to 3no can be reached according to 
model calculations [15]. Unfortunately, in the case of heavy-ion collisions such measurements 
face a number of problems related to the expansion dynamics and the spatial inhomogeneity 
of the produced dense matter. Due to the collective expansion, the baryon density n depends 
on time so that obviously a time average of the mass shift is to be expected, i.e. the mass 
shift effect is distributed over an interval related to the temporal change of n. Also spatial 
variations of n will make the situation more complicated. 

To extract the wanted information on the in-medium meson spectral function in heavy- 
ion collisions one has also to take into account that the di-electron spectrum dN/dM 2 as 
a function of the di-electron invariant mass M is actually a convolution of the considered 
spectral function and the probability to create the vector meson with given energy and mo- 
mentum. The latter one usually drops down very fast with the energy making the resulting 
di-electron spectrum very steep. This can make the identification of the possible mass shift, 
superimposed to collision broadening and further in-medium broadening, in particular for 
the p meson, rather difficult. An additional problem arises for the u meson because of its 
comparatively long life time which, however, may change in a dense medium [11]. To mea- 
sure an in-medium mass shift of the u resonance one should obviously separate lepton pairs 
produced by decays in the dense medium from such ones which originate from decays after 
the disintegration of the dense system. 

The aim of the present note is to explore the extension of the usual analysis of the 
invariant mass spectrum dN/dM 2 to the more informative double-differential spectrum 

dN/dM 2 dM]_, where M± = \J M 2 + Q\ is the transverse mass of the pair with transverse 

— * 

two-momentum Q± perpendicular to the beam axis. As shown in [16] within a schematic 
fireball model, such an extension with suitable kinematical cuts allows to avoid the major- 
ity of the above mentioned difficulties in extracting information on the in-medium vector 
meson spectral function. Here we present results of our calculations of double-differential di- 
electron spectra in heavy-ion collisions at SIS18 energies obtained within a transport model 
of Boltzmann-Uehling-Uhlenbeck (BUU) type. 



II. DI-ELECTRON PRODUCTION AND VECTOR MESON SPECTRAL 

FUNCTION 

As elementary processes for the e + e~ production we employ pion annihilation 7r + 7r~ — > 
p — > e + e~, direct vector meson decays V — > e + e~ with V = p, u>, pion - nucleon and nucleon - 
nucleon bremsstrahlung, and Dalitz decays of mesons, A and N*. Pions are produced solely 
in baryonic resonance decays, and these resonances are created either in baryon - baryon 
collisions via BB — > BB with B — N, A, N* or in meson - baryon collisions MB — > B 
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with M = n,p,a, while the "direct vector mesons" emerge in the elementary reactions 
BB — > VNN or 7ri? — > ViV. Our treatment of these processes is described in some detail 
in [18]. To be specific, let us list our perturbative treatment of the p and u channels. 

We use the parameterization [17] of the total cross section <r vac for 7r + 7r~ — > p — > e + e~ 
in vacuum 



o- vac (M) = ^ r^y ,/i-M ^ (i) 

3 VM7 V M 2 (M 2 — m 2 ) 2 + (rhpTp) 2 V ' 

with m p = 775 MeV, rh p = 761 MeV, T p = 118 MeV to generate the in-medium spectrum 
by 

a(M) = o- vac (M) — p \ - / V ; -, 2 

V ; V ; m 2 (n = 0)A(M;m p (n = 0),r p (n = 0))' KJ 

where the spectral function is 

1 m r 

A(M; m p , T p ) = -— p p r ^ . (3) 

p p/ 7T (M 2 - m 2 ) 2 + (m p r p ) 2 

with mass dependent width 

r>) = rro (B) 1-5^23=1 . (4) 



\Jrn 2 p (ri) - Ami 



Below we employ certain parameterizations of T p °\n) and m p (n). 

Similar to the case of 7r + 7r~ annihilation, the distribution of in-medium decays of directly 
produced p mesons into e + e~ follows in local density and pole mass approximation from 1 

^f- = a p (Vs) 2M A(M; m p (n), r») Br(p - e+e") (5) 

with the branching ratio of Br(p — > e + e~) = F(p — > e + e~)/T(p, tot), so that all in-medium 
effects are related to the modification of the vector-meson spectral function. We take the 
default value of Br as in vacuum, i.e., assume that the e.m. decay width and the total width 
change in the same manner. Such a choice of the branching ratio allows us to avoid an 
artificial change of the relative contributions of p and u> mesons (see below for the case of 
the density dependent u> mass shift). In Eq. (5), cr p (y/s) is a parameterization of the total 
cross section. For NN — > NNp we use a p — F P S > where x = (y/s — m p (n))/GeV and 
S = 1, F p = 1.5 mb, £ = 1.4, y — 1, z = 2. In case of NA (AA) collisions the scaling factor 
S = 0.75 (0.5) is employed. The nN — > iVp reaction is similarly parameterized, but F p = 1.5 



1 This contribution is of sub-leading order (see below), therefore, an estimate of the maximum 
yield at pole position is sufficient here. 
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mb, f = 0.018, y = 2.2, z = 3.5. The scaling factors S here are 0.375 for 7rA++[A-], 0.125 
for 7T+A , 7T-A+, and 0.25 for tt°A. 

The e + e~ decay spectrum of u mesons follows from 



f = ^ dtexp{-T UJ (n) 1 t}2MA(M;m UJ (n),T UJ (n))T LU (n)Br(uj^e + e-), (6) 



where the path through matter with local density n(t, x) is followed. 7 is the Lorentz factor 
according to the velocity, and Br(uo — > e + e~) denotes the vacuum branching ratio. The 
production cross section a w is given by the above parameterization with F w = 0.36 mb, 
£ = 1.25, y = 1.4, z = 2 for the reaction NN — > NNu with the same scaling factors as for 
the irN — > Np recation and corresponding isospin channels. The reaction ttN — > Noo follows 
from = 1.38 mb, £ = 0.0011, y = 1.6, 2; = 1.7, where the appropriate scaling factors 5 
are 0.5 for 7r°iV[iV*], 1 for 7r + n or n~p, and all others as above. 



According to our experience [18], the di-electron spectra roughly scale when changing 
the mass numbers of the colliding nuclei. We choose, therefore, the scaled cross section 
da sc&l /dM = irR 2 A- 2 dN(b = 0)/dM with R = A 1/s 1.124 fm. Since we are going to 
demonstrate the advantage of certain observables for the case of strong medium effects, we 
focus on central collisions of gold nuclei at a beam energy of 1 AGeV. 

First we use the parameterization T p °\n) = aT p , m p (n) = m p — Am p with constant 
values of a and Am p . This is schematic and aimed at verifying the previous finding [16] 
that a strong in-medium broadening of the vector meson widths smears out the p meson 
peak in the invariant mass spectrum, while the double-differential spectrum still allows to 
identify the p peak. This set should also illustrate the separate effects of mass shifts. While 
in [16] a fireball model with thermal and chemical equilibrium was used to calculate the di- 
electron spectra, here both these conditions are not longer required. Moreover, to contrast 
the present study with the one in [16] we switch off any rescattering of the p and uj mesons 
once created thus not allowing thermalization and chemical equilibration of these species. 

In Fig. 1 we exhibit the resulting double-differential cross section da scal / dM dM± for the 
parameters a = 3 and Am p = 0. Let us first consider the invariant mass spectrum arising by 
the full Mj_ integration (see lower panel in Fig. 1). We are interested in the invariant mass 
region M > 400 MeV, where obviously the pion annihilation channel n + ir~ — > p — > e + e~ 
dominates over Dalitz decay and bremsstrahlung contributions. The oj peak sticks out, but 
the chosen broadening is so strong that any peak structure, attributable to the p meson by 
pion annihilation, disappears. Instead, a knee is reminiscent of the p. This is the general 
expectation that a spectroscopy of the p meson becomes impossible due to broadening [6]. 
However, selecting an appropriate M± window, say 850 MeV < M±_ < 950 MeV, the p 
peak becomes clearly visible, even for the assumed extreme broadening (see upper panel 
in Fig. 1). Due to phase space weighting, the p peak position appears slightly shifted to 
smaller mass. At the high-M wing the direct p contribution and the pion annihilation 
yield compete. The M± window around 900 MeV seems most suitable, as it represents a 
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compromise of suppressing unwanted background, not cutting off kinematically the p peak, 
and still allowing for sufficient statistics. 

Next let us consider the effect of a mass shift, which is strongest for the p meson according 
to many predictions. Fig. 2 demonstrates for a = 3, Am p = 150 MeV that the above found 
features are persistent also in case of quite drastic mass shifts. The upper panel in Fig. 2 
illustrates, indeed, that such a strong mass shift of the p meson is clearly visible even in 
case of the strong broadening in the spectrum da scal /dM g?M_i_|m ± =850---950McV, while for the 
chosen parameters the invariant mass spectrum does not exhibit any p peak, due to the 
assumed strong broadening (see lower panel in Fig. 2). 

Figs. 1 and 2 show that for identifying the uj meson still the invariant mass spectrum 
is enough: For in-medium changes covered by the parameterization of set (i) there is no 
obvious advantage by a particular M± selection. The important point here is that the uj 
meson rides on the pion annihilation background. 

Note that in Figs. 1 and 2 the distributions are according to the spectral function Eq. (3) 
and the phase space; in line with the model assumptions (no density dependence of Amy 
and Ty, no rescattering) , there is no need of an explicit propagation of the spectral function. 

Let us now turn to temporal and spatial variations of the density which cause an ad- 
ditional distribution of the vector meson decay strength. Also here, the usefulness of the 
double-differential cross section with appropriate kinematical cuts still applies. We focus on 
an illustration of the uj mass shift. Since the uj in vacuum is a narrow resonance, often it is 
assumed to survive the dense stages during a heavy-ion collision and to decay after freeze- 
out with its vacuum parameters. Therefore, the uj in-medium spectroscopy is considered as 
more useful in reactions of elementary projectiles with a nucleus under special kinematical 
conditions in which the uj is roughly at rest and decays inside the static matter. We would 
like to analyze here the prospects of the uj in-medium spectroscopy in heavy-ion collisions. 

In focusing on the uj meson, we leave for the moment being the p meson parameters 
unattached (i.e., take vacuum parameters) and choose m w (n) = — dm^n/no, r u (n) = 
T w (n = 0) + ST UJ n/no with Sm^ = 70 MeV and ST^ = 50 MeV entering the spectral function 
of the uj which is of the form of Eq. (3). Given the still sufficiently narrow width of the 
in-medium uj we propagate the peak according to local conditions and distribute it at each 
time instant with the spectral function Eq. (3) and sum up the decay products, cf. Eq. (6). 
These procedure is acceptable only for sufficiently narrow resonances (see discussion in [19]). 
Since the number of uj mesons is normalized, an increase of the width means a reduction 
of the peak height. For the selected parameters, the uj peak for the above in- medium 
parameters has a contribution from in-medium decays and vacuum decays. The first one 
shows clearly the down shift and the smearing due to the density dependence, see lower panel 
in Fig. 3. The vacuum peak sticks not longer so much above the pion annihilation yield. 
The in-medium and vacuum uj contributions compete with the p contribution from 7r + 7r~ 
annihilation. To increase the relative contribution of in-medium uj decays, in [20] it has been 
proposed to select \ow-Q± di-electrons which stem from very slow uj mesons which decay 
in the medium, in contrast to fast uj mesons which traverse quickly the medium and decay 
outside. This feature is also seen in our simulations for heavy-ion collisions in the upper 
panel in Fig. 3, where we selected di-electrons with Q± < 100 MeV: the vacuum contribution 
becomes significantly less, while the in-medium uj contribution is clearly seen as a shifted 
bump riding on the 7r + 7r~ annihilation background. This is in contrast to the invariant 
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mass spectrum which is integrated over all Q± (see lower panel in Fig. 3) where the strong 
competition between the wanted uj meson signal and the 7r + 7r~ annihilation background 
makes a conclusion on a in- medium uj mass shift difficult. Therefore, the measurement of 
the double-differential di-electron spectrum at low values of Q± offers much better chances 
to verify a uj meson mass shift. Further cuts, e.g. setting an additional narrow repidity 
window around mid-rapidity, suppress too strongly the yield but do not enhance efficiently 
enough the in-medium contribution. 

With respect to the latest HADES e + e~ measurement in carbon - carbon collisions, we 
mention that in our simulations of central collisions C + C the Q± dependence is similar 
to the case of Au + Au. However, for C + C the "in-medium u shoulder" appears below 
the 7r + 7r~ channel, since the in-medium uj contribution is attenuated due to the shorter life 
time and the smaller volume of compressed nuclear matter. The vacuum decay contribution 
dominates, even for small values of Q±. 

Finally, we would like to comment on the above choice of Br = Br vac , which assumes 
that T(cj — > e + e~) and r(u;,tot) scale in the same manner in medium. In contrast, often 
an in-medium change of r(a;,tot) is assumed, and T(uj — > e + e~) is kept constant (cf. [15] 
for a discussion). In the latter case the in- medium contribution is noticeably suppressed, as 
evidenced in Fig. 4. Clearly, this issue, as the treatment of dynamical in-medium effects of 
p and uj on a common footing, needs further investigations. 

IV. SUMMARY 

Our results can be summarized as follows: (i) The spectrum da/dM dM± for 850 MeV 
< M ± < 950 MeV as a function of M is dominated in the region 400 MeV < M < 800 
MeV by the channel 7r + 7r~ — > p — > e + e _ with a pronounced p bump. This gives a good 
chance to measure the in-medium p meson spectral function by a double-differential di- 
electron spectrum in heavy-ion collisions at SIS 18 energies, (ii) By selecting di-electrons 
with small enough transverse momentum Q± < 100 MeV one can isolate the in- medium 
uj decay contribution which signals the in-medium mass shift and a broadening of the uj 
meson. This Q± dependence may be useful for future uj spectroscopy, (iii) Bremsstrahlung 
and Dalitz decays are small in the region M > 400 MeV, for the considered beam energies. 

Our findings, which are not to be understood as quantitative predictions, also show that 
fairly independently of the specific form of the vector meson spectral function the double- 
differential di-electron spectrum with suitable cuts allows a better access to the in-medium 
effects and, therefore, can be an appropriate measurement strategy in HADES experiments. 
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FIG. 1. Contributions to the di-electron spectrum da scal /dM for various bins of M± for central 
collisions Au(l AGeV) + Au. A resonance broadening by a factor a = 3 is assumed as in [16] for 
both p and uj\ no mass shift. 
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FIG. 2. As in Fig. 1 but with additional mass shift Am p = 150 MeV. 



9 



0.001 



I 

> 

CD 

o 



03 

o 
w 


T3 



I 

> 
CD 



05 

o 
w 

D 
T3 



0.0001 



1e-05 



0.01 



0.001 



0.0001 



1e-05 




0.5 0.6 0.7 0.8 0.9 1 
M [GeV] 




0.5 0.6 0.7 0.8 0.9 1 
M [GeV] 



CO 



tot 
w vac 
^Vned 
7l + K - p 

p deacys 



FIG. 3. Contributions to the di-electron spectrum da scal / dM for various bins of Q± for central 
collisions Au(l AGeV) + Au. The parameters are ST^ = 50 MeV and bm^ = 70 MeV. For an 
orientation, also the p contributions are shown, but with vacuum parameters. 



10 



0.001 



I 

> 

CD 

o 



03 

o 
w 


T3 



I 

> 
CD 



05 

o 
w 

D 
T3 



0.0001 



1e-05 



0.01 



0.001 



0.0001 



1e-05 



Q T < 100 MeV 



_i i_ 




0.5 0.6 0.7 0.8 0.9 1 
M [GeV] 




0.5 0.6 0.7 0.8 0.9 1 
M [GeV] 



CO 



tot 



-vac 
+ ^med 
71 71 - p 

p decays 



FIG. 4. As in Fig. 3, but with an to branching ratio as described in text 
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